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Abstract. Formulating a quasiclassical approach we determine the cross section for 
the complete four-body break-up of the lithium ground state following single photon 
absorption from threshold up to 220 eV excess energy. In addition, we develop a new 
classification scheme for three-electron ionizing trajectories in terms of electron-electron 
collisions, thereby identifying two main ionization paths which the three electrons 
in the ground state of lithium follow to escape to the continuum. The dominant 
escape paths manifest themselves in a characteristic "T-shape" break-up pattern of the 
three electrons which implies observable structures in the electronic angular correlation 
probability. This break-up pattern prevails for excess energies so low that the Wannier 
threshold law a oc E a describes already the triple ionization cross section, whose 
predicted value a = 2.16 we can confirm quantitatively. 
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1. Introduction 

The broad interest in multi-electron photoionization processes is due to the fundamental 
role they play for understanding electron correlation induced by the long-range Coulomb 
forces. Identifying the different paths the electrons follow to escape to the continuum 
is essential in uncovering a variety of fundamental phenomena, from pattern in highly 
differential cross sections over interference phenomena to the energy dependence of the 
total ionization cross section. The theoretical treatment of multiple ionization processes 
is highly complex with no analytic solution. In the energy domain, the difficulty is that 
one has to account for the correlated motion of the electrons in the asymptotic form 
of the final continuum state. In the time domain, this difficulty can be avoided at the 
expense of propagating the fully coupled few-body Coulomb problem in time. 

For three electrons and sufficiently high photon energies this has been achieved in [1] 
and for even higher energies one may use approximate schemes [2] . Experimentally, the 
first data on triple ionization appeared only recently [3,4], compared to data on photo 
double ionization which date back to the late sixties of the last century [5] . 

In the current work, we present a theoretical study of the total triple ionization 
cross section of lithium over a wide range in energy, from threshold up to 220 eV 
excess energy, and predict characteristic features in the electronic angular correlation 
probability. While the total triple ionization cross section has already been measured [3] 
and compares favorably with our results, the electronic angular correlation probability 
is not known experimentally. However, it should be measurable with state of the art 
experimental techniques. 

Given the obstacles in the theoretical description stated above we are only able 
to achieve these results by formulating the four-body break-up process quasiclassically. 
This implies classical propagation of the Coulomb four-body problem using the classical 
trajectory Monte Carlo (CTMC) phase space method. CTMC has often been used 
to describe break-up processes induced by particle impact [6-9] with implementations 
differing usually in the way the phase space distribution of the initial state is constructed. 
We use a Wigner transform of the initial quantum wave function for the initial state, 
and this is why we call our approach "quasi" -classical. Naturally, the electron-electron 
interaction is treated to all orders in the propagation, and any difficulties with electron 
correlation in the final state are absent, since the method is explicitly time-dependent. 

While the classical results follow from a large numerical effort, they still allow 
for a detailed analysis of the trajectories in terms of their physical properties. We 
will demonstrate that the triply photoionizing trajectories can be organized in groups 
according to the respective sequence of electron-electron collisions. From the emerging 
scheme we identify two main paths that lead to triple ionization from the Li ground 
state. The group to which an ionizing trajectory belongs is identified in an automated 
process which warrants a transfer of the classification scheme to more than three-electron 
atoms without technical difficulties. Physically, the nature of the collision scheme also 
promises a generalization to more electrons. 
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Furthermore, being quasiclassical, our approach naturally addresses the energy- 
regime close to threshold. This energy regime has been traditionally of particular interest 
since the slow electron escape allows for large interaction times resulting in pronounced 
interactions among the escaping electrons. We could confirm the Wannier threshold 
law a oc E a for the four-body break-up of Lithium by single photon absorption with 
a = 2.16 in the energy range of 0.1-2 eV as detailed in [10]. 

The paper is organized as follows: In section 2 we explain the theoretical approach, 
in section 3 we present the most important results, the triple ionization cross section 
and the electronic angular correlation probability. Section 4 introduces the classification 
scheme of ionizing trajectories which can explain the electronic angular correlation 
probability in terms of the dominant "T-shaped" pattern of the three escaping electrons. 
Section 5 is dedicated to a final discussion and a summary. 

2. Quasiclassical theory of photoionization 

We formulate the triple photoionization process from the Li ground state (ls 2 2s) as 
a two step process [2,11,12]. First, one electron absorbs the photon (photo-electron). 
Then, due to the electronic correlations, redistribution of the energy takes place resulting 
in three electrons escaping to the continuum. We express the above two step process as 

a 3+ = ^absP 3+ , (1) 

where <7 a b s is the total absorption cross section and P 3+ is the probability for triple 
ionization. In what follows, we evaluate P 3+ and use the experimental data of Wehlitz 
[13] for (T a bs- Equally well, we could use a theoretically calculated (T abs [14] which is 
easy to obtain following the approach of [15]. To compute P 3+ , we first assume that 
the photo-electron is a ls-electron. It absorbs the photon at the nucleus (r x = 0), 
an approximation that becomes exact in the limit of high photon energy [16] and 
is in analogy to the successful description of photo double ionization in two-electron 
atoms [12]. The photon could also be absorbed by the Li 2s-electron. However, the 
cross section for photon absorption from a Is orbital is much larger than from a 2s 
orbital as investigated in [17] for photoionization of an excited He(ls2s) atom. Hence, 
we can safely assume that the photo-electron is a Is electron which significantly reduces 
the initial phase space to be sampled. Also, by virtue of their different character the 
electrons become practically distinguishable and allow us to neglect antisymmetrization 
of the initial state. We denote the photo-electron by 1, the other Is electron by 2 and 
the 2s electron by 3. Immediately after photon absorption, we model the initial phase 
space distribution of the remaining two electrons, Is and 2s, by the Wigner transform 
of the corresponding initial wavefunction ip(r 1 = 0,r 2 ,r 3 ), where are the electron 
vectors starting at the nucleus. We approximate the initial wavefunction as a simple 
product of hydrogenic orbitals ^(r,) with effective charges Z i} to facilitate the Wigner 
transformation. The Zi are chosen to reproduce the known ionization potentials Jj, 
namely for the 2s electron Z 3 = 1.259 (I 3 = 0.198 a.u.) and for the Is electron Z 2 = 2.358 
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(J 2 = 2.780 a. u.). (We use atomic units throughout the paper if not stated otherwise.) 
The excess energy, E, is given by E = E^ — / with E u the photon energy and / = 7.478 
a.u. the Li triple ionization threshold energy. For a given E, the Wigner distribution W 
has an energy spread and it is only its expectation value that is equal to E [18]. Since 
near E — energy conservation is vital we enforce it by restricting the Wigner functions 
for the individual electron orbitals to their respective energy shell [19]. Following these 
considerations, the initial phase space distribution is given by 

p (V) = ^V5(r 1 )5(e 1 + h - u) J[ (r u Pi )8(ei + h) (2) 

i=2,3 

with normalization constant jY . 

We determine the triple ionization probability P 3+ formally through 

P 3+ = lim / dr p3+ exp((t - t abs )^f cl )p(r), (3) 

u r abs 

where the classical Liouvillian Jzf c i is defined by the Poisson bracket {H, } with the full 
Coulomb four-body Hamiltonian H [20] and the propagation begins at the time t^s of 
photoabsorption. The projector V 3+ indicates that we integrate only over those parts of 
phase space that lead to triple ionization. Equation (JHJ) amounts to propagating electron 
trajectories using the classical equations of motion (CTMC). Regularized coordinates 
[21] are used to avoid problems with electron trajectories starting at the nucleus. We 
evaluate P 3+ by weighting each triply ionized trajectory by the initial phase space 
distribution and adding the contributions [12]. 

3. Experimentally accessible observables: Total triple photoionization and 
angular correlation probability 

Observables calculated classically or semiclassically tend to better approximate the 
"exact" value if they represent a quantity averaged over as many degrees of freedom as 
possible. The reason for this is that classically forbidden mechanisms such as tunneling 
and pronounced interference effects are less likely to play a prominent role if the 
dynamics is averaged over many degrees of freedom. Secondly, for classical calculations 
based on Monte Carlo methods, one "counts" events, similarly as in the experiment. 
This means, only (randomly sampled) trajectories, whose final phase space values fall 
into a certain bin, contribute to the observable. As in the experiment, statistics of the 
contribution is an important factor and therefore, more integrated observables are easier 
to determine compared to highly differential ones. Interestingly the classical domain of 
validity is complementary to that of approximate quantum calculations, e.g., with the 
so called 6C wavefunction [22] where the fully differential cross section is calculated and 
more integrated observables require numerically expensive integrations. 

For the present case, this implies that our focus is the triple ionization cross 
section (which is the "most integrated" three-electron observable) with our statistics 
also allowing the evaluation of single differential cross sections. Here, we present for 
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reasons which will become clear in section 3 an unusual but observable quantity, the 
electronic angular correlation probability, which answers the question: How likely is it 
in triple ionization to find two electrons under a certain angle 91 

3.1. Triple photoionization cross section 

We have already described in section 2 how to determine the triple photo ionization 
cross section o~ 3+ which is shown in figure ^i. Note that no fitting parameters are used 
to obtain a 3+ . However, one may object that we use the total photo cross section 
extracted from another calculation or the experiment for the curve in figure Hence 
we present in figure \I]p the triple ionization probability P 3+ which we calculate directly 
and which can be also obtained directly from experimental data using the relation 

p3+ = jR 3+ /(1 + jR 2+ + jR 3+ ); (4) 

where R 3+ and R 2+ are the experimental triple-to-single and double-to-single 
photoionization ratios, respectively [3,4,13]. Figure also illustrates the numerical 
effort involved in computing P 3+ . For example, to obtain 10 3 triple photoionizing 
trajectories at E = 0.9 eV with P 3+ ~ 10~ 7 one has to evolve 10 10 trajectories with the 
CTMC method [6]. 

Considering the approximations we had to make to handle the four-body problem, 
the agreement with the experimental data is remarkably good, starting near threshold 
where we can confirm the classically expected behavior of the cross section [10] according 
to Wannier's theory [23,24]. The agreement extends beyond the maximum of the cross 
section where our results also agree with the data points recently obtained in an ab- 
initio calculation [1]. Hence, our classical approach with an approximate initial quantum 
wavefunction apparently captures the relevant correlations among the three electrons. 
For very high excess energies (currently not considered) the triple photoionization 
can not be adequately described by our quasiclassical formulation. However, one can 
describe the process using the Born approximation. 

3.2. Angular correlation probability 

As already indicated, with only 10 3 triple ionization events out of 10 10 trajectories, at 
most a single differential observable can be determined. Interesting with respect to 
dynamical correlation among the three electrons and experimentally accessible is the 
angular correlation probability 

C{d) =lim V / dr P a + exp((t - t ab8 )^ cl )p{Y)8{e ij {t) - 9) (5) 

i>j=l J tabs 

with 

dij{t) = arccoslpiitfpj (t) / (pi^pj^t))] (6) 

which depends only on the relative angle 0y between any pair of ionized electrons in 
the three electron escape. Formally, this is easily achieved within a CTMC approach. 
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Figure 1. (a) Triple photoionization cross section obtained by multiplying the triple 
photoionization probability from the present calculation (o) with the total photo cross 
section from [f 3] in comparison to the experiment [3,4] (experimental results (•) shown 
with error bars), (b) shows triple photoionization probability P 3+ as a function of 
excess energy: present calculation (o); experiment (•) [3,4,13]. 



However, the number of events is very small and we need to bin the observable over 
10° respectively, as shown in figure 121 We see for the higher excess energy (E = 220.5 
eV) a broad distribution with a maximum near 90°. This might have been expected 
from impulsive binary collisions of the fast photo electron (peak at 9 = 90°) where the 
shift to slightly larger 9 indicates the (small) influence of the Coulomb repulsion. At 
E = 0.9 eV, the situation is drastically different: A double-hump structure emerges with 
peaks at 90° and 180°^ This can be interpreted as "T-shaped" structure which the three 
outgoing electrons form, where two of them leave along a line towards opposite sides 
and the third one leaves perpendicularly to the line. The origin of this double-hump 
structure and its interpretation as a T-shape configuration of the escaping electrons will 
become clear after the analysis of the electron collision sequences in the next section. 

4. Classification of triple ionization dynamics in terms of an 
electron-electron collision scheme 

The definition of electron collisions in a multi- electron system 

The (dipole) coupling to the electromagnetic field is a single-electron operator. Hence, 
one electron absorbs the photon, which gets thereby annihilated. Before this photo 
electron has left the atom it must transfer part of its energy to the other electrons to be 
ionized within a very short time. How this happens, directly or indirectly, is a question 
which is difficult to ask quantum mechanically, but is a natural question in classical 
mechanics, where the electrons undergo soft collisions mediated by Coulomb forces. In 

\ Convolution with the volume element sin# leads to an appearance of the 180°-peak at a smaller 
angle, see section |4~31 
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Figure 2. The angular correlation probability C(9) /P 3+ (E) denned in JSJ for E = 0.9 
eV (•) and for E = 220.5 eV (o), normalized to the respective triple ionization 
probability (P 3+ (E = 0.9eV) = 7.89 x 10~ 8 , P 3+ (E = 220.5eV) = 6.43 x 10~ 5 ). 
The number of events N has been binned over 10° at 8j — 5° + jl0°, j — 0,1, ... , 17. 
For the solid line, see section l4~3l 



contrast to billiard balls they can indeed transfer energy among themselves if the nucleus 
is within the reach of the Coulomb potential to absorb the recoil momentum. Moreover, 
in a two-electron atom, it is most likely that a single collision among the two electrons 
occurs which does not provide a lot of insight. However, in three-electron atoms, the 
situation is far more complex and a priori the nature of the collisions is not clear: Is 
triple ionization mediated mostly by a single collision involving all three electrons or 
does it happen sequentially with a sequence of momentum transferring two-electron 
collisions? If the latter is the case, is there a pattern of preferred sequences? Does a 
characteristic collision pattern of classical trajectories lead to observable consequences? 

Before we can answer these questions we first have to define what we call 
a momentum transferring electron-electron collision along a trajectory with time 
dependent electron positions r-j(t), i = 1,2,3. The term responsible for momentum 
transfer between electrons i and j is their Coulomb repulsion V^r^) = , = — rj. 
Hence, we identify a collision between electron i and j (ij) through the momentum 
transfer 



under the condition that V(rjj(tfc)), k = 1, 2 are local minima in time with t<i > t±, while 
Tij = \ri — Tj\. This automatically ensures that the integral of (JJJ) includes the "collision" 
with a local maximum of V{vij(tk)) at a time t\ < tM < ^2- During the time interval 
t\ < t < t 2 , all four particles interact with each other. Hence, the definition (JJJ) is only 
meaningful if the collision redistributes energy dominantly within the subsystem given 
by the two-electron Li + -Hamiltonian Hjj of the nucleus and the electrons % and j involved 
in the actual collision. This is indeed the case since the energy in the subsystem H^- 
changes little over the collision, i.e., dHij/dtdt <C E, where E is the total energy of 
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Figure 3. Energies in time for three-body subsystems for a triple ionizing trajectory 
with (12,13) electron-electron collision sequence at E = 220.5 eV. The maxima of 
V(ri2) and V(ri3) occur at times 0.017 a.u. (0.41 attoseconds) and 2.22 a.u. (54 
attoseconds), respectively. In (a) the solid line is the energy of the H12 subsystem, 
while the dashed line is the interaction V(ri2). In (b) the solid line is the energy of 
the H13 subsystem, while the dashed line is V(rxa). The insets are for short times t. 



the Li-Hamiltonian. We illustrate the latter statement in figure fusing as an example a 
triple ionizing trajectory labeled by the sequence of (12,13) electron-electron collisions. 
During the 12 collision (see figure Ek), V(r 12 ) undergoes a sharp change while H 12 
changes smoothly. Thus, during the 12 collision, the V(ri 2 ) potential energy is primarily 
redistributed in the Li + subsystem of electrons 1 and 2. Similarly, during the 13 collision 
(see figure Eb), V(ri 3 ) undergoes a sharp change while if 13 changes smoothly. Thus, 
the V(ris) potential energy is primarily redistributed in the Li + subsystem of electrons 
1 and 3. On the other hand H13 and H12 undergo a sharp change during the 12 and 13 
collisions, respectively. This should be the case, since during these times it is the energy 
of the H 12 and H 13 subsystem, respectively, that is conserved. Note that the higher 
the excess energy is, the more impulsive the electron-electron collisions are, since the 
collision time becomes shorter and shorter compared to the time the bound electrons 
need to orbit around the nucleus. 

4-2. Dominant collision sequences 

For the majority of triple ionizing trajectories, we register at least two electron- 
electron collisions. For automated identification of the collisions sequences we need 
a sensitivity threshold to register only the important collisions for the triple ionizing 
trajectories. This is done individually for each trajectory by forming the maximum 
D = maXj^dDjjl} and normalizing each collision according to Dij = |Dy|/D. A 
collision is only registered if > 5. The resulting classification shown in figure 0] does 
not sensitively depend on the exact value of 5 which we have chosen to be 5 = 1/8. 

According to our classification scheme the photo-electron transfers energy to the 
other two electrons through two main sequences of collisions in about 84% of all triple 
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Figure 4. Probability relative to all triple ionizing trajectories classified according to 
the sequence Si of electron-electron collision ij, see text; si = (12,13) (o), s 2 = (12,23) 
(.), s 3 = (13,12) (□), S4 = (12,13,23) (A), s 5 = (23,12,13) (a). 



ionization events. In the first main pathway to triple ionization, the si=(12,13) sequence, 
the photo electron (1) knocks out, successively, electrons 2 and 3. In the second main 
pathway, the $2= (12,23) sequence, the photo-electron 1 first knocks out electron 2. 
Then, electron 2 knocks out electron 3. It is easy to understand how the relative 
probability of these two processes is changing as a function of excess energy. The process 
Si=(12,13) has the highest weight for low excess energy where the photo-electron (after 
photon absorption) is still slow enough to easily transfer energy, first to electron 2 and 
then to electron 3. For higher energy, the competing process S2=(12,23) takes over 
because the photo-electron is still so fast after its first collision with electron 2 that 
the interaction with the more loosely bound 2s electron 3 is small. Rather, it is more 
probable that electron 2 transfers part of the energy it has gained after the initial 12 
collision to electron 3 through the 23 collision. 

Generally speaking, the more collisions a sequence has, the less is its weight which 
can be seen in figure 0] comparing the weight of S4 and S5 with s\ and s%. However, 
there are also exceptions, as illustrated by S3 = (13, 12). That this sequence does not 
carry large weight is obvious since it is not very likely that the photo-electron collides 
first with (2s) electron 3 and later with the more tightly bound (Is) electron 2. 

4-3. The influence of the collision scheme on the angular correlation probability 

One would expect that the collision scheme with its two dominant pathways to the 
continuum, the sequences s\ and S2, has a prominent influence on collisional observables. 
Indeed, as we will show, it explains the T-shaped configuration (see figure EJ) of the three 
outgoing electrons which manifests itself in the angular correlation displayed in figure El 
To this end, we have investigated the angular correlation C{9) in more detail by asking 
which angle 9ij the individual electron pairs ij (i.e., 12, 23, and 13) form upon leaving the 
nucleus. This allows us to determine the characteristic C (6) a specific collision sequence 
Si produces. We show the result in figure El for the two dominant sequences Si and S2- 
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Figure 5. Sketch of the T-shape structure when the three electrons escape to the 
continuum through the si process (a) and through the S2 process (b). 



One sees that the electron pair which undergoes the last collision in a sequence (13 for 
Si and 23 for s 2 ) leaves towards opposite sides with an angle of 9 = 180°. On the other 
hand, the pair which collides first (12 for both sequences) forms an angle of 9 = 90°, 
and so do the pairs which do not undergo an explicit collision (23 for s± and 13 for s 2 . 
The solid and dashed lines in figure El are fits with the functions 

Cf{9) = cYsnA0sin0 (8) 
CjJ(0) = cj'sin /3 n^/2sin^, (9) 

respectively, where the c* J and (3 are fitting parameters, while sin 9 comes from the line 
element of integration, sm9d9. Taking the direction of the electron leaving in the plane 
perpendicular to the pair escaping back to back as a reference, it is easy to see that the 
average width in the distribution about 180° should be twice the one about 90°, i.e., 
~ 2(3± = 2(3. This is indeed the case and all curves in figure El are fitted well by 
(3 = 16.8. The fits to the individual curves in figure El together produce the curve shown 
in figure El according to 

C{9) = sin9[c ± su/(fl) + c,, sin 2/3 (#/2)] (10) 

with 



c± = ci 2 ( Sl ) + ci 3 ( Sl ) + ci 2 ( S2 ) + ci 3 ( S2 ) 
q = c f( Sl )+cf(s 2 ). 



(11) 
(12) 



For the two dominant collision sequences s\ and s 2 the T-shape geometry implies that 
the inter-electronic angle 90° is produced twice (by the first pair of electrons colliding 
and by the pair which does not collide) while only the pair involved in the second and 
last collision produces a peak at 180°. Since the peaks at 90° and 180° for the angular 
correlation of individual electron pairs have roughly the same height (see figure EJ) this 
should lead in the full angular correlation probability to a peak at 90° which is twice as 
high as the one near 180°. As one can see from figure [2] this is indeed the case. 

The T-shape is the consequence of the sequential process of two collisions which 
dominates triple ionization since both main pathways Si and s 2 contain two collisions. 
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Figure 6. Angular correlation probability for E — 0.9 eV as in figure^ but between 
two specific electrons 12 (open square/solid line), electrons 13 (filled circle/dotted 
line) and electrons 23 (open circle/dashed line). Part (a) is for trajectories from the 
Si sequence, part (b) for those from the S2 sequence, for details, see text. 



Widely unnoticed the T-shape configuration was already mentioned in [25] for small 
excess energy in a mo del- calculation of electron impact ionization of Helium. The reason 
that this T-shape like escape has not received any attention since then lies probably in 
the fact that its connection to the fundamental organization of the triple escape in 
characteristic collision sequences of the electrons was not known until now. 

4-4- Evolution of the angular correlation distribution with excess energy 

As already indicated in figure 121 the T-shape is lost towards higher energies where 
impulsive collisions dominate and the electrons are so fast that 9 = 180° is not reached. 
On the other hand, approaching E — > we expect the symmetric triangular escape 
according to Wannier. In our present calculation we see a clear tendency for the 
transition to the Wannier configuration only at the lowest excess energy (0.1 eV) we are 
able to calculate, see figure Hf. This could mean that, at least concerning differential 
observables, the Wannier regime in triple ionization has an extremely short range in 
energy. It could, however, also mean that our approach, treating the photo electron 
very differently from the other electrons, favors asymmetric break up of the electrons 
and therefore a very low transition energy between symmetric and asymmetric electron 
configurations. 
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Figure 7. Angular correlation probability C{9) for excess energies of E =10 eV (a), 
6.2 eV (b), 2 cV (c), 0.9 eV (d), 0.5 eV (e) and 0.1 eV (f). C{9) has been binned over 
20° at 10°, 30°, . . . , 150°, 170° to reduce the statistical error. Its value for the total 
triple ionization probability does not exceed 3% in (a-e) and 10% in (f). 
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5. Discussion and conclusions 

To summarize, we find classically that in triple ionization of lithium by single photon 
absorption the three electrons escape to the continuum through a dominant T-shape 
configuration for excess energies that are in the range of validity of the Wannier law 
(a = 2.16) [10]. This T-shape configuration gives rise to a double hump structure in 
the inter-electronic angular distribution of the escaping electrons at 90° and 180°. We 
have explained this surprising double hump structure in terms of a novel classification 
scheme that is built upon momentum transferring electron-electron collisions. 

In the framework of many body perturbation theory (MPT) different ionization 
processes have been known for many years [26] . However, as the name already says, the 
perturbative character of this approach makes it applicable at high excess energies only 
with the electron-electron interaction treated to first or second order. The classification 
scheme we have developed emerges from the full classical dynamics which includes the 
electron-electron interaction to all orders. Surprisingly, our scheme shows that we can 
describe the triple photoionization process as a sequence of electron-electron collisions 
for energies close to threshold, where the Wannier theory becomes valid. 

The sequence of collisions relevant for lithium involves only three-body Helium-like 
subunits (nucleus and two electrons) at one instant of time. Since Coulomb systems 
interact via two-body forces only, it may well be that this scheme holds for more than 
three electrons in an atom. At the same time this would imply that not a two-electron 
atom but a three-electron atom is the fundamental system whose understanding allows 
one to access multi-electron ionization dynamics in the future. Our scheme might also 
guide the way to a quantum mechanical analysis along the lines of [27]. 

Furthermore, our study may offer valuable insight into the triple photoionization 
process in connection to double ionization by electron impact. For double ionization, the 
relationship between electron impact ionization of He + and a quasiclassical formulation 
of the double photoionization from the He ground state has already been established 
[11, 12]. Both processes differ only slightly, namely, in the energy scale set by the 
respective bound electron. We believe that our quasiclassical formulation of the triple 
photoionization of the Li ground state is a process very similar to double ionization 
of the excited states ls2s 1 ' 3 S' of Li + by electron impact. The target state involved 
is the ls2s Li + , since in the triple photoionization process the photo-electron is a Is 
electron knocking out the remaining two, Is and 2s, electrons. Experimental/theoretical 
results on double ionization of Li + by electron impact are needed to establish whether or 
not the two three-electron escape processes are indeed similar. Such results would also 
elucidate how the spin symmetry — not accounted for in our classical calculation — affects 
the double ionization by electron impact from the excited states ls2s 1,3 S of Li + . 

An experimental investigation of the double hump structure we predict for the 
inter-electronic angular distribution for small excess energies is very desirable. It could 
confirm the existence of the T-shaped escape of the three electrons and would thereby 
support our classification scheme of collisions as well as the validity and limitations of 
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a classical approach to the four-body Coulomb problem. 
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